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ABSTRACT 

The membership of Matsunaga's variable 1, a carbon-rich, mass-losing, Mira variable, 
in the globular cluster Lynga 7 is discussed on the basis of radial velocities. We con- 
clude that it is a member, the first known C-Mira in a globular cluster. Since such a 
variable is expected to have an age of ~ 1 — 2 Gyr and an initial mass of ~ 1.5 solar 
masses, we conclude that this star must be the product of a stellar merger. 
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1 INTRODUCTION 

In the course of an extensive infrared survey of Galactic 
globular clusters for variable stars, Matsunaga (2006) dis- 
covered a long period Mira variable near the centre of the 
globular cluster Lynga 7. The variable has a period of 551 
days, a large infrared amplitude (AK — 1.22 mag) and very 
red colours (mean (J—K) = 4.1 mag) indicative of mass-loss 
and a circumstellar dust shell. Sloan et al. (2010) obtained 
SPITZER mid-infrared spectra of the variable showing it 
. to be a carbon star with strong circumstellar dust emis- 
sion from SiC and MgS particles as well as absorption from 
gaseous acetylene. They estimate a mass loss of 2.5 x 10~ 7 so- 
lar masses per year. A few carbon rich objects are known in 
globular clusters, but these have generally been interpreted 
as cluster examples of the CH stars in the Galactic halo field. 
These latter stars are believed to be mass exchange binaries. 
However no carbon Mira has previously been detected in a 
globular cluster and Matsunaga's variable, VI, raises a num- 
ber of problems related both to stellar evolution and to the 
nature of Mira variables. 

In the present paper we derive the radial velocity of 
VI and discuss its cluster membership. We then discuss the 
nature of VI and its relation to the cluster population. 



2 OBSERVATIONS 

A medium dispersion (resolution 2.4A) spectrum of VI 
Lynga 7 was obtained on 2012 May 11 (HJD 2456058.6) 
with the RSS spectrograph on the SALT telescope at SAAO, 



* Based on observations made with the Southern African Large 
Telescope SALT. 
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Sutherland. The spectral region covered was 6000 to 7100A. 
Spectra of the planetary nebula Hen 2-146 and the carbon 
Mira V650 CrA (P = 332 days) were also obtained with the 
same equipment and settings. 

Fig. 1 shows the spectrum of VI and that of V650 CrA. 
The spectra match each other well. In both cases the spec- 
trum is dominated by molecules of carbon compounds (par- 
ticularly CN) and is typical of a carbon-rich Mira (see for 
instance the atlas of Barnbaum, Stone & Keenan (1996)). In 
addition VI shows strong Ha emission. This is characteristic 
of C-Miras at certain phases. 

Using the comparison arc we obtain for the planetary 
nebula Hen 2-146 a radial velocity^] of 56±4kms~ 1 in good 
agreement with the catalogue value of 61 ± 4km s _1 (Du- 
rand at al. 1998) . Measured against the comparison arc the 
Ha emission in VI gives 17km s -1 with an estimated un- 
certainty of ~ 2km s^ 1 . The night sky lines in this spec- 
trum yield 6±2kms~ 1 . Thus we adopt the 11km s _1 as our 
best estimate of the velocity from the Ha emission. This 
emission originates in shock waves in rising material in the 
stellar atmosphere. It therefore has, predominately, a nega- 
tive velocity with respect to the true radial velocity of the 
star. The best measure (if available) of a C-Mira velocity 
is obtained from CO lines in the millimetre region, which 
arise well above the photosphere (see, e.g., Barnbaum 1992 
and references therein). Data discussed in Menzies et al. 
(2006) show that on average Vel(CO(mm)) = Vel(Emission) 
+ 16km s" 1 . Thus our best estimate of the stellar velocity 
is 27km s^ 1 . The uncertainty in this figure is not easily de- 
termined, but is likely to be 5km s _1 or more. 

A radial velocity was obtained for VI by cross- 



Throughout we use heliocentric radial velocities. 
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Figure 1. Comparison of the Lynga7-Vl spectrum (upper panel) 
with that of the reference carbon star, V650 CrA (lower panel). 
The Ha line in Lynga7-Vl has a peak count of 7880 e — (offscalc 
at relative count of 1.64), while the peak count at 7020A in V650 
CrA is 1.3xl0 5 e - . The counts go to zero in the gaps between the 
CCD chips. 



correlating its spectrum with that of V650 CrA, after mask- 
ing out the strong Ha emission line. A radial velocity of 
S^kms- 1 was adopted for V650 CrA (Walker 1979). Sev- 
eral measures were made using all or parts of the spectral 
range and we adopt a value of 40 ± 10km s _1 as the best 
mean estimate. In Menzies et al (2006) an average correc- 
tion of — 4kms _1 of C-Mira absorption velocities to the CO 
(mm) standard system was adopted, thus giving a estimate 
of 36 ± 10km s -1 for the stellar velocity. The absorption line 
velocity of a C-Mira varies in a cycle by ~ 13km s~ (San- 
ford 1944). This will affect both VI and V650 CrA and the 
cross correlation velocity is thus likely to be somewhat less 
certain than that derived from Ha emission. 



3 CLUSTER MEMBERSHIP OF VI 

The position of VI given by Matsunaga (2006) 
(16 h ll m 02 s .O - 55°19'l4"(J2000)) together with the 
cluster centre given by Goldsbury et al. (2010) leads to a 
distance of the variable from the centre of A(RA) = —13.5 
and A(Dec) = —10 arcsecs well within the core of the 
cluster (core radius 54 arcsecs, half-light radius 72 arcsecs 
(Harris 1996)). 

Radial velocities of 9 cluster stars were measured by 
Saviane et al (2012). They find a mean heliocentric veloc- 
ity of 22±3kms _1 . Earlier Tavarez & Friel (1995) obtained 
6 ± 15km s -1 from 4 stars. The high standard error in this 
case being due to an uncertainty in night to night zero 
points. In section 2 we obtained a radial velocity for VI from 
Ha emission 27± ~ 5km s -1 with a somewhat less certain 
value of 36 ± lOkrns -1 from a match to V650 CrA. Given 



the uncertainties we consider the agreement with the cluster 
velocity to be satisfactory. 

In addition to this, the velocity of VI differs from that 
expected for a field C-rich Mira of its period. Long period 
C-Miras in the solar vicinity are disc objects showing evi- 
dence of differential Galactic rotation. The velocity disper- 
sions are moderate, 22±4kms _1 in the Galactic radial direc- 
tion and 18 ± 3km s -1 in the direction of Galactic rotation, 
at the relevant period. (Feast et al. 2006). These disper- 
sion include observational scatter and may be slight over- 
estimates. Assuming a flat rotation curve (O = 220km s -1 
and Ro = 8.0kpc and with the local solar motion from Feast 
& Whitelock (1997), the predicted radial velocity for VI at 
the distance predicted by a bolometric period-luminosity re- 
lation (see section 4) is — 81km s -1 , quite different from the 
value observed. Also for a member of a disc population with 
a relatively small velocity dispersion VI is quite far from the 
Galactic plane at its predicted distance (445 pc). 



4 DISCUSSION 

Saviane et al.(2012) determined a metallicity for Lynga 7 of 
[Fe/H] = -0.57 ± 0.15 on the scale of Carretta (2009) and 
Sarajedini et al. (2007) show that an HST colour- magnitude 
diagram is well fitted by that of 47 Tuc ([Fe/H] =-0.76). 
Lynga 7 has been claimed to be among the oldest globulars 
(Marin- Franch et al. 2009). Sarajedini et al. estimate (m — 
M)„ = 14.55 and E(B-V) = 0.78 for Lynga 7. This is based 
on a fit to the 47 Tuc c-m diagram and adopting 13.40 and 
0.055 for the latter. There have been some differences in the 
distance estimates for 47 Tuc (see Gratton et al. 2003) and 
the uncertainty in the modulus of Lynga 7 is probably not 
less than ~ O.lmag. 

For VI, Matsunaga obtained mean magnitudes of J = 
11.35, H = 9.08 and K s = 7.25 (J — K s = 4.10) on the 
IRSF system which is close to the 2MASS system (Kato et 
al. 2007). Converted to the SAAO system using the trans- 
formations given by Carpenter (2001 and web update) and 
applying a reddening equivalent to E(B — V) = 0.78 on the 
Cardelli system (Cardelli et al. 1989) we obtain rriboi = 9.60 
using the BC K -{J-K) relation of Whitelock et al. (2006) . 
A main uncertainty is probably the conversion from IRSF 
(assumed same as 2MASS) to SAAO. Thus with the Sara- 
jedini distance modulus the bolometric absolute magnitude 
of VI is —5.0. The value predicted from the the period- 
luminosity (Mboi) relation of Whitelock et al.(2009) is —5.2. 
This relation is based on LMC C-Miras with bolometric cor- 
rections calculated in the same way as for VI and adopting 
a modulus of 18.5 for the LMC. Taking into account the 
various uncertainties including the intrinsic scatter in the 
period-luminosity relation, the observed and predicted val- 
ues are in satisfactory agreement. Note that VI is 0.9 mag 
fainter at K than predicted by a period-luminosity relation 
in K (Whitelock et al. 2008), confirming that the star is 
surrounded by an obscuring dust shell. 

Galactic kinematics (Feast et al. 2006) and, particularly, 
the presence of mass-losing carbon Miras with periods ~ 500 
days in intermediate age (1-2 Gyr) clusters in the Magellanic 
Clouds (e.g. Nishida et al. 2000) strongly suggests that the 
initial mass of a C-Mira such as VI was about 1.5 solar 
masses. VI also matches properties of these cluster variables 
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well in absolute magnitude and colours. On the other hand, 
Padova isochrones (Marigo et al. 2008, Giradi et al. 2010 and 
website) for 10-13Gyr globular clusters indicated that the 
initial mass of their evolved stars is about 0.8 solar masses. 

Since the discovery of a carbon-rich star (a "CH star" ) 
in u) Cen by Harding (1962) several such stars have been 
discovered. These are taken to be examples of mass exchange 
binaries of which many are known in the general halo field. 
In such cases the amount of mass exchanged is likely to 
have been relatively limited and sufficient only to pollute 
the stellar atmosphere. Sloan et al. (2010) suggest that VI 
could be a CH type object. However, with a predicted initial 
mass of ~ 1.5 solar masses this is not possible. The formation 
of a star of this mass in Lynga 7 requires the merger of two 
stars near (or beyond) the main sequence turn of point of 
the cluster. 

Sharina et al. (2012) have found a carbon star mem- 
ber in the globular cluster NGC 6426. It has quite strong 
molecular bands but is relatively blue in JHK. Nothing is 
known about its variability. The authors discuss whether it 
is a mass exchange binary or a merger and are not able to 
come to a certain conclusion (since the mass is not known). 
A somewhat similar carbon rich star was found in the glob- 
ular cluster NGC6405 (M15) by Cote et al. (1997), which 
they consider to be a typical CH star, i.e. a binary with a 
white dwarf secondary which earlier donated mass to the 
observed star. 

The nature of blue stragglers in the c-m of globular clus- 
ters has been much discussed since they were first identified 
by Sandage (1953). It now seems that there are two types of 
blue stragglers, mass exchange binaries and stellar mergers 
(see Ferraro et al. 2009 and references therein). An initial 
discussion of the evolution of those blue stragglers which are 
merged stars has been given by Sills et al. (2009). Further 
theoretical work might allow one to determine whether the 
properties of VI require that the merger components be near 
the main sequence or whether mergers of giant branch stars 
(whose internal chemical composition has been changed by 
evolution, e.g. helium enrichment) would also result in a 
typical long period C-Mira. It would be interesting to know 
whether there are blue stragglers in Lynga 7. However, radial 
velocities or proper motions are required for this (Milone et 
al. 2012). In any case objects like VI must be very rare since 
their life time as mass-losing AGB variables is quite short . 



South Africa (NRF). The observations reported in this paper 
were obtained with the Southern African Large Telescope, 
SALT. 
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5 CONCLUSIONS 

The radial velocity of Matsunaga's variable VI is consis- 
tent with membership of the globular cluster Lynga 7. Its 
absolute bolometric magnitude and colours are typical for 
a long period mass-losing carbon-rich Mira variable. Since 
such variables have initial masses of ~ 1.5 solar masses and 
ages of ~ 1.5 — 2.0Gyr, we propose that VI is the evolu- 
tionary product of a stellar merger in the cluster ~ 1.5Gyr 
ago. 
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